LARGE DEVIATIONS AND LINEAR STATISTICS FOR POTENTIAL 
THEORETIC ENSEMBLES ASSOCIATED WITH REGULAR CLOSED SETS 
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o 

3 ■ Abstract. A two-dimensional statistical model of N charged particles interacting via logarith- 

mic repulsion in the presence of an oppositely charged regular closed region K whose charge 
£f) ' density is determined by its equilibrium potential at an inverse temperature fi is investigated. 

When the charge on the region, s, is greater than TV, the particles accumulate in a neighborhood 
of the boundary of K, and form a point process in the complex plane. We describe the weak* 
1 limits of the joint intensities of this point process and show that it is exponentially likely to find 

^ ^ 1 the process in a neighborhood of the equilibrium measure for K. 

" 1. Introduction 

In two-dimensional electrostatics, charged particles are identified with points in the extended 
complex plane. The potential energy of a system of two like charged particles located at z,w G C 
is proportional to — log \z — w\. More generally, if zi, z%, . . . , zn are the locations of N identically 
charged particles, then {z\, . . . , zm} determines the state of the system and the potential energy 
^ , of this state is given by 

2^ ; - lo si z « ~ z m\- 

The energy is minimized when the particles are all at oo. In order for the system to be found in a 

^ * . state where the particles are at finite positions, there needs to be a potential (or other obstructions) 

which repels the particles from oo. We represent this field by V so that the interaction energy 
CS| . between a particle located at z and the field is given by V(z). The total potential energy of the 

system comprised of the N particles in the field is given by 

N 

E(zi, . . . , z N ) = ^2 V(z n ) - ^2 log \z n - z m \. 

n—l m<n 

The system is assumed to be in contact with a heat reservoir so that the energy of the system is 
variable, but the temperature is fixed. In this setting, j3 denotes the reciprocal of the temperature, 
and the Boltzmann factor for the state {z\, . . . , zjy} is given by 



X 



N 



71=1 



This quantity gives the relative density of states, so that the probability (density) of finding the 
system in state {z±, . . . , z?f} is given by 

Z- N \ v e-^--^\ Z NAV =[ e-e E ^>->^dA® N ( Zl ,...,z N ). 

Jc N 
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Let u := Yl n =i be the empirical measure associated with a state {z\, . . . , z/v}. In this work 
we take V to be minus the equilibrium potential for some regular closed region K and study the 
following two questions. First: what is probability to find the system (empirical measure uj) close 
to a given Borel measure on the complex plane (large deviation principle); second: what is the 
limiting behavior of the marginal probabilities of the probability density function for this system 
(linear statistics). 

These questions have been obviously considered before. In [17, 18, 2], see also [10, 1], the large 
deviation principle was shown for the case of the external field V being identically +oo of the real 
line (thus, all the charges are confined to the real line) and satisfying limi-ri^oo (log \x\ — eV(x)) = 
— oo for any e > 0. This work was further extended to the case V(z) = \z\ 2 without confinement to 
the real line in [3, 8]. In [9], the large deviation principle was shown to holds for particles restricted 
to the unit circle, i.e., V is continuous on the unit circle and V = +oo of it. 

As to the linear statistics, the case of particles confined to the real line and V being polynomial 
of even degree with positive leading coefficient was studied in [11]. The case where particles are 
restricted to a compact subset of the complex plane interacting in the presence of a continuous 
field was treated in [6]. No confinement case with regular fields satisfying V(z) > (1 + e) log(l + |z|) 
for some e > was considered in [7]. 

2. Main Results 

2.1. Potential Theoretic Setting. For any probability Borel measure on C, say v, set 



:= f log^^^ 2 (z, U ) 
J \z-u\ 



to be its logarithmic energy (negative free entropy), where &v® n (zi, , . . . , z n ) := dv(z\) ■ ■ ■ &v(z n ). 
For any compact set K the logarithmic capacity of K is defined by 

cp(K) := exp < — inf I[u] 

I supp(^)CX 

It is known that either cp(-ftT) = (K is polar) or there exists the unique measure wk, the 
logarithmic equilibrium distribution on K, that realizes the infimum. That is, 



(1) cp(if) = exp | - I[u K ] 

The logarithmic potential of ujk , that is, 

V""(z) := / log-^—dw K (u), 



is equal to I[uik] quasi everywhere (up to a polar set) on K and is at most as large everywhere in the 
complex plane. The set K is called regular with respect to the Dirichlet problem if V UJK = I[cok] 
everywhere on K. The Green function with a pole at infinity for the unbounded component of K c 7 
the complement of K, is defined by 

g K :=I[w]-V*. 

It is a non-negative harmonic function in K° \ {oo} with a logarithmic singularity at infinity. 
Moreover, it is zero q.e. on K and is, in fact, continuous if K is regular. 

Let s and TV be two parameters such that s > N . We always assume that s and N scale in 
such a fashion that the limit of s~ 1 N as N tends to infinity exists and define the following energy 
functional 

(2) ItW\ := I[v] + ^[ g K dis, I := lim s^N, 
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where it is understood that Iq[v] = oo for all v such that supp(i/) n K c ^ 0. 

Proposition 1. Let K be a compact set with connected complement which is regular with respect 
to the Dirichlet problem. For all £ £ [0, 1] , it holds that 

(3) I[u K ] < h[v] 

for any compactly supported probability Borel measure v ^ ljk- 

Clearly, h[v\ = I [is] for any measure v supported on K as gx = on K. Hence, Ii[lok] = I\wk] 
and therefore ujk is the unique minimizer of the weighted energy functional Ii for any i £ [0, 1]. 

2.2. A Model for Random Configurations. Let if be a compact set with connected comple- 
ment which is regular with respect to the Dirichlet problem. In this paper we investigate random 
configurations whose joint density is given by 



(4) n N , s ,fs(z 1 ,...,z N ) := — cxpi-(3sy^g K (z n )} \ \z r , 



N 1 

-/j*x>(*0f n i 

n— 1 ) jn<n 



where (N, s, j3) is a triple of numbers such that 

(5) /3(s- N + l) > 2 + c 

for some fixed cq > 0, and Zn >s ,p is a normalizing factor that turns Oj\r,s,/3 into a probability 
density function. Clearly, 

(6) z n, s ,i3= / exp{-Psy^g K (z n )\ fdA 9N (z!,... ,zn) 



\ exp<^ -(3s Vff K (z n ) \ 



where dA stands for the Lebesgue measure on C. 

Besides connections to electrostatics and random matrix theory, the results below are also 
motivated by number theory. Let K be such that cp(K) = 1 and p be a polynomial. The Mahler 
of p with respect to K is defined by 

M K {p) := exp i / log \p\dw K I = cxp i a p ^ g K (z) >, 

where a p is the leading coefficient of p. When K = D, it is known that dwjf(z) = ^r|dz|, and 
therefore M := Mp is simply the classical Mahler measure. In [4], the bound for the number of 
polynomials with integer coefficients of degree at most TV such that M(f) < const, was derived. 
The main term of the asymptotics for this bound came from Zm,n +1,2 defined by (6) with K = D. 
Moreover, it was shown that Z^^^ is a rational function of s with poles at every positive integer 
less or equal to N. The cases where K is an ellipse and E = [—2, 2] were also investigate in [14, 15]. 

2.3. Large Deviation Principle. Let r\ = {771, . . . , tjn} be a random configuration chosen ac- 
cording to the law Ojv.s,^- That is, the probability that rjj <S O, j £ {1,...,JV}, is equal to 
J QN flN,s,0dA® N for any open set O C C. To any such configuration we associate the empirical 
measure defined as 

N 



N ^ J ^ 

fe=i 

where S z is the classical Dirac delta distribution with the unit mass at z. 

Let v and /i be two probability Borel measures on C. Then the distance between them is defined 

by 



dist(i/, /i) = sup 
/ 
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where the supremum is taken over all functions / that are bounded by 1 in modulus and satisfy the 
Lipschitz condition with constant 1 on supp(i^) Usup(/x). For measures supported on a compact 
set it holds that dist^, v n ) — > as n — > oo if and only if v n A v as n — > oo, where A stands for 
the convergence in the weak* topology of measures (J f&v n — > J fdv for any continuous function 

/)• 

Theorem 2. Let K be a compact set with connected complement which is regular with respect to 
the Dirichlet problem and such that K = K°. Given (5), it holds that 

(7) Jim^logZ w = -|/K]. 
Moreover, for any probability Borel measure v, supp(z^) C C, it is true that 

(8) lim lim log Prob {dist (u, u n ) < e} = - - ( I e [v] - I[u K ] 

e— N— >oo 7V 2 \ 

Equations (3) and (8) yield that the probability to find utr, in a small enough neighborhood of v 
is subexponentially small if Ie[v] = oo and is exponentially small if Ii[v\ < oo and dist (i/,uk) > 0. 

The asymptotics in (7) can be improved if we restrict the attention to Jordan domains with 
smooth boundary and /3 = 2. 

Proposition 3. Let K be a Jordan domain whose boundary dK is a Jordan curve of class 1 C 1,ot , 
a > 1/2. Then 

(9) logZ w , s , 2 = -N(N + l)I[u K ] + 9{l)N + O(logTV), 

where 0(x) = log7r + 1 + — x) log(l — x), which is a continuous increasing function on [0, 1] 

with values 9(0) = log7r and 9(1) = logir + 1. When s = oo, the term 0(logN) can be replaced by 

0(1). 

2.4. Linear Statistics. The ?i-th marginal probability of ftjy s p, n <G {1, . . . , N — 1}, is defined 

by 

(10) ^\jz u ...,z n ) := / n N ^(z 1 ,...,z N )dA^ N - n \z n+1 ,...,z N ). 

J C N-n 

Let r} be a random configuration chosen according to the law £In,s,p and oj v be the corresponding 
empirical measure. Then uj^ can be considered as a point process on C. It is known that the 
joint intensities of this point process are equal to (N\/(N — n.)!)Qjy s (3 . That is, ifOi,...,O n are 
mutually disjoint subsets of C, then 



~[w v (O k ) 



.fe=i 



AH 



(N-n)\J 0lX ... x0n 



where E[-] denotes the expected value of a random variable. The following theorem describes the 
weak* behavior of the measures oo. 

Theorem 4. Under the conditions of Theorem 2, it holds that 

lim f filftlpdA*" = f ffop 

for each f € C&(C"), n G N, where Cb(C") is the Banach space of bounded continuous functions 
on C n . 



^The arclcngth function of dK is continuously diffcrcntiablc as a periodic function on the real line and its derivative 
is a-H61dcr continuous. 
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Theorem 4 in particular implies that 



lim E 

jV->oo 



fdcj 



= / fdoJK, f G C b (C), 



where r] N = {fyf , . . . , t/n} 1S a random configuration chosen according to the law £In,s,/3 since 



e [/«)]= //nS!.^ 



for any j e {1, . . . , TV}. More generally, let {ji, . . . ,j n } C {1, . . . , N} be a set of distinct indices 
and fx,..., f n G C 6 (C). Set . . . , z„) := JlLi /*(**)• Thei1 / e C b (C") and it holds that 

where the inner sum is taken over all possible combinations of n— k+ 1 coordinates being equal. Ob- 
serve that the first integral on the right-hand side of the equality above converges to J^)! =1 J /fcdwjf 
as N — > co by Theorem 4. This observation yields the following corollary to Theorem 4. 

Corollary 5. Let f G C&(C) and k, m G N. Under conditions of Theorem 2, it holds that 



lim E 

AT-KX) 



3. Proofs 



/dwif 



Proof of Proposition 1. Let i/ be such that dist(i/, wjr) > 0. Without loss of generality we may 
assume that v has finite energy. Recall that gx = on K. Thus, if supp(V) C K, then h[v\ = 
> I[ujk] by the very definition of luk- Otherwise, consider 

where v\K" is the balayage of v\ K c onto K relative to K c , [13, Sec. II. 4]. Then it follows from [13, 
Thm. II. 4. 7] that v has finite energy as well and V v < V u + J gudv in K c . Integrating both sides 
of this inequality against v and using Fubini-Tonclli's theorem for the left-hand side, we get that 

J V v du < I[v] + j g K dv. 

In fact, it also true that V v = V u + J gxdv everywhere on K as K is regular with respect to the 
Dirichlct problem, [13, Sec. II. 4]. Therefore, 

< I\v\ + 2 J g K dv = h[v] < I t [u]. 

The desired conclusion now follows from the fact that supp(P) C K and therefore < I\y\- □ 

Put, for brevity, Wk '■= er 9K and define 

Sn(w K ) ■= SUp < ] [ \z n - Z m \ JJ^k" 1 ^") 1 ■ ■ -i z n) G C W > . 

\,m<n n ) 



(ii) 

Lemma 6. Let {X\ 



, Xn} be any configuration satisfying 



Sn(w k ) = J [ | An - A TO | JJl 



Z- 1 



(A„), 



then {Ai, . . . , A at} c K . 
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Proof. Set 

Pk, N {z) :=w^-- 1 (z)l[(z-X n ). 

Then ||P fc) iv||c = |-Pfc,Jv(A fc )| for each k G {1, . . ., N} since 

8n{w k ) = ||-Pfc,w||c | An - A m | TJ Wx{\ n ). 

m<n, n.m^k n^k 

As deg(Pfe i Ar) = N — 1, % = 1 on K, and wk(z) < 1 for z <E AT C , Bernstein- Walsh inequality [13, 
Thm. III.2.1] yields that 

\P k . N (z)\ < \\Pk, N \\K < \\Pk,N\\c z e 

Hence. P^jv attains its maximum on K and therefore Xk G if for all k G {1, . . . , A}. □ 
Let P be a compact set satisfying all the conditions of Theorem 2. Define 



(12) E m := \z G E\ dist(z,<9P) > — \, m G N. 

I m J 

Lemma 7. Let E be as described and is be a probability Borel measure supported on E with finite 
energy. Then there exists a sequence of probability Borel measures {v m } such that supp(f m ) C E m , 
dist {y m , v) — > and I[v m ] — > I [is] as m — > oo. Moreover, if v = lue, then v m can be taken u>E m . 

Proof. For each m the measure v can be written as 

v = V \E m + v \E°\E m + V\ dE . 

By the very definition of the sets E m it holds that E m —\ C E m and E° = U m E m . Hence n m (P° \ 
E m ) = and therefore |f|.E°\.E m | — > as m — > oo. If v is supported only on the boundary of P, 
v = v\qe, we set v m := v m , where v m is the balayagc of v\qe onto E m relative to E^. Otherwise, 
v\E m is not a zero measure for all m large enough and therefore we can define for such m 

V m '■= Ot m V\E m + Vm, 

where a m := 1 + \v\E°\E m \l\ v \E m \ is chosen so \v m \ = 1, a rn — > 1 as m — > oo. Clearly, supp(^ m ) C 
E m - 

The sequence {v m } has a weak* limit point, say v*. Since v*{B) = v(B) for any compact set 
B C E° , it holds that v* EO = v\ E ° by the interior regularity of Borel measures. Therefore, 

v* = v ]E o + v* dB , 

where Iz^^l = |^|asl ~ \v m \ and vT gE is a weak* limit point of {pm\- Let us show that v* gE — v\qe- 
This will imply that v* = v and therefore v m — > v as v* is an arbitrary weak* limit point of {u m }. 

Let A C N be a subsequence such that v m — > v* SE as m — > oo, m G A. Since supp(^*g B ) C dE, 
it holds that 

V v ' aE (z)= hm y p ">0), zgP°. 
On the other hand, by the very properties of balayage [13, Thm. II.4.7], it holds that 

V Vm {z) = V v \°»{z)+c m , z G E° m , 

where c m = J gE m <lv\dE an d 9E m is the Green function for P^ with pole at infinity. Since 
E m+ i C E^, the maximum principle for harmonic functions applied to gE m ~~ 9E m+1 yields that 
the functions gE m form a decreasing sequence on dE. Thus, c := lim m _ >00 c m is well-defined. 
Hence, we have that 

V"\ob = yv\8E + c on E ^ 
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Since logarithmic potentials are continuous in the fine topology by the very definition of the latter 
[13, Sec. 1.5], the equality above can be extended to every fine limit point of E° . As E is regular 
with respect to the Dirichlet problem in E c , the fine closure of E° coincides with E [13, Thm. A. 2.1]. 
Thus, 

V V * BE = V VvdE + c on ED( supp(v l9E ) U supp(i/f aB )) . 
As both measures have finite energy and therefore are C-absolutely continuous, it holds that 
v \dE = v \dE an d c = by [13, Thm. II. 4.6]. 

To show that I[u m ] — > I[v] as to — > oo, it is enough to prove that limsup m ^ 00 I[v m ] < I[v] as 
the opposite inequality follows from the principle of descent [13, Thm. 1.6.8]. It is a straightforward 
computation to get 

V = V m - (a m - l)l/|jgo + a. m V\ E °\E m + V\l)E - v m , 

where, again, the two middle terms containing a m are not present when v = v\oe- Thus, since 
V v \be _ yVn, > _ Cm in £ [ 13) Thm n.4.7], it holds that 

V > V v ™ - (a m - 1)V^° + a m V v \^\^ - Cm 

> V Vm - (a m - 1)1/^° - a m \v ]E o\ Em \ logdiam^) - c rn 

on E, where we also used an observation that V a + \a\ logdiam(iJ) > on E for any positive 
measure a supported on E. Recall that c m , |i / |s o \E m | — ^ and a m — > 1 (whenever present) as 
to — > oo. Hence, integrating both sides of the last inequality against v and taking the limit superior 
of the right-hand side yields 

I[v] >limsup / y I/ " 1 d^ = limsup / V h dv m . 

m— ¥oo J m—*oo J 

Using the estimate from below for V v once more with the same caveat concerning a m , we get that 



> limsup fl[v m ] - (a m - 1) f V U \ E ° dt 



The desired inequality follows now from the fact that the integrals J V V \ E ° &v m are uniformly 
bounded above. Indeed, notice that 

V v \*° < V" + - \v ]B o\) logdiam(£;) 

and therefore we get by what precedes that 

limsup / ~Y V \ B ° du m < limsup / V v du m + (\v\ - \v\eo\) logdiam(^) 

< I[v] + - |i/| S o|)logdiam(E). 

Finally assume that v = loe- Since supp(ws) C dE, it holds that v m is simply the balayage of 
uj e onto E m relative to Ef n . Then V Um = V^ E + c m = I[uj e ] + c m q.c. on E m [13, Thm. II.4.7]. 
However, the latter property uniquely characterizes equilibrium measures [13, Thm. 1.3.3]. □ 

Let v be a compactly supported probability Borcl measure. Define 
(13) Av e := a 6 dA, a e (z) := — / l s (u - z)dv(u), 



where l e is the indicator function of the disk {\u\ < e}. So defined, v e is also a probability measure 
with compact support. 

Lemma 8. It holds that dist (u, v e ) < e and I[v e ] — > I[v] as e — >• 0. 
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Proof. Let / be a Lipschitz continuous function on supp(V e ) with Lipschitz constant 1 for some 
e > 0. Observe that supp(V) c supp(V e ). Then 

fdu e -Jfdu = J (/(*) - f(u)) le ^~ Z) dv(u)dA(z) <J\z- u\ le ^~ z) dv(u)dA(z) < e 

by the Fubini-Tonelli theorem and since \f(z) — f{u)\ < \z — u\ < e for z, u € supp(V e ), \z — u\ < e. 
Hence, dist (i/, z^ e ) < £ and therefore z^ e A ^ as £ — > 0. 
If I[v\ =oo, then 



/mM :=_/min{M,log ] ^}d I 



diverges to infinity as M — > oo by the monotone convergence theorem. Since the functional Im 
is defined with respect to a continuous kernel, it follows from the weak* convergence of measures 
that 

liminf I[u E ] > lim ImWb] = hiM- 

Because the inequality above is true for any M, the limit of I[v e ] as £ — > diverges to infinity as 
well. The rest of the lemma is the content of [7, Sec. 3.2]. □ 

For further use, let us state the following trivial modification of the principle of descent [13, 
Thm. 1.6.8] for empirical measures. 

Lemma 9. Let {^rj N } be a sequence of empirical measures that converges weak* to some probability 
Borel measure v . Then 

I\v] < liminf I.[u VM ], /.[«„„]=/ log j— — dw®*(«, u). 

Proof. Since to rlN — > u>, it holds that w® 2 — > v® 2 . Hence, it follows from the monotone convergence 
theorem and the continuity of min |m, log \ z l _ u \ j on C 2 that 

IM = lim / min \ M, log , 1 , \&v® 2 {z,u) 

M->oo J [ \Z — U\ J 

= lim lim / min i M, log r 1 du;® 2 (z, «) 

M-Hx> JV-S-oo J \ b \z-ll\j nNX ' 

< lim liminf ( IJu)„J + — ) = liminf I Jw„J. 

Af->oo jV->oo AT y N^oo N 

□ 

Lemma 10. Let v e be as in (13). Then there exist configurations rj N ' £ such that mmj^k \n^' £ — 

Vk' e \ — c'/v^Vj dist [y € . oj^n.e) = O (A^" 1 / 4 ), and I* [uj^n.c] — > I[v e ] as N — > oo, for some constant 
c' that depends on i/ e but does not depend on N . 

Proof. For convenience, set M := [A^ 1 / 2 ] . Since v e is absolutely continuous with respect to d^4, 
there exist real numbers x\ < xi < ■ ■ ■ < xm and a positive constant such that the vertical 
strips 

M 



S 3 : : : x 3 < Re(z) <x j + -^\, je {1,...,M}, 



are mutually disjoint and supp(^ e ) C ] J^L 1 Sj. Clearly, &m is bounded above by diam(supp(^ e )) 
and below by the 1-dimcnsional Lebcsgue measure of the projection of supp(i/ e ) on the £-axis. 
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Moreover, it holds by (13) that 

f c M diam(supp(i/ e )) 
u e ( Sj ) = ^ du £ < — . 

In particular, the inequality above implies that 

Mj := [u e (Sj)N\ ^dN 1 / 2 

(all the constants Ci appearing from now on depend only v e ). Furthermore, as v e is a unit measure, 
it holds that 

M 

N-M <Y^ M j < N - 
i=i 

Let now be the greatest ordinate such that v e {Sj PI {z : > Im(z)}) = 0. Then for each 
k G {1, . . . , Mj} there exists yj.k+i such that 

v e (Rj,k) = N^ 1 , R jtk ■= Sj n {z : y jtk < Im(z) < yj,k+i} ■ 
Observe that by (13) it holds that 

7T£ 2 M ire 2 1 c 2 

Set 

rj jtk := Xj + iy jtk , k G {1, . . . , M 3 + 1} and j G {1, . . . , M}. 

Then the collection r\ := {r)j t k} contains at least N and at most N + M points and by the very 
construction 

min{6M,c 2 } d 



min \r] j k - r^j > —rj- > 

for some d > that depends only on v e as the constants bu are uniformly bounded above and 
below. 

Let R := \jf =l U*^ R j>k . Then 

W=[EE / + 2 EE / J logr-^d,f 2 (,,.)=:^I, + 2^/,, i . 

To estimate the first sum notice that the monotonicity of the logarithm and the choice of the 
rectangles Rj >k yield 

Mj Mi 

i > 1 VVi 1 

Moreover, by the very choice of the points rjj jk , we have that \r)j,k — > c 2\k — llN^ 1 / 2 . Thus, 
as the width of each R^ k is bja/M, we deduce that 



I > fc, 



diam (Rj jk U Rjj) < 



Hence, using the fact that 



Mi Mi . Mj Mj 



1 x ^ \ ^ / C3 \ 1 x - \ ■> C3 C4 Mj 

2iV2^^ og l + a + 1 - fc) 2 ) - ^m 2 ^ ^ a + 1 - k) 2 - ~T W 

k=l l=k v v ' / fc=l i=fc v ; 
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and that exactly the same estimate holds for I <E {1, . . . , k — 1}, we get that 

1 Mj+1 Mj+1 1 1 Mj 1 M 

and respectively that 

M 1 M Mj+1 

To estimate the second sum, fix j,i € {1, . . . , M} such that j < i. It can be readily observed that 
for each fee {1, . . . , Mj} there exists 1%. E {1, . . . , Mj} such that 

\V],k+i ~ Vi.i + ^>m/M\, I < l k , 

\Vj,k - Vi,l+1 + bAl/M\, I > l k . 



diam (Rj t k U Ri,i) 
Therefore, it holds that 



1 A / \Vj,k+i - I <h, 

diam(i? Jifc U R u ) < ( 1 



I \Vj,k - vu+i\, i>h- 

By the above inequality and since Mj < ciN 1 / 2 , we deduce that 

1 Mi ( lk ~ l 1 Mi 1 \ 1 c 2 

It can be readily verified that the sequence {lk}^f=i is non-decreasing and therefore the pairs of 
indices in the double sums above never repeat themselves. Thus, 

for some sequence {(k m J m )}. So, using the fact that \r]j t k m ~ Wi,l m \ — c'N' 1 / 2 , we get that 

„ x - x - 1 log TV 

Combing the estimates for Jj and ^ J^j, we derive that 

I[v e ] = I[u e]R }+ f V v ' +v '^dv e]R . > limsup (h[ur,\ - c 7 ^|^ + / V^+^«dv E \ Rc ) = lim sup 7, 

where the last equality follows from the fact that v e \ R c < N~ x / 2 and since the potentials V Ue+u '\ r 
are uniformly bounded on supp(^ £ ) as the measures v e -\-v B \ R are absolutely continuous with respect 
to dA with uniformly bounded densities. 
To show that 

I[ v e] < liminf h[u) v ], 

N^oo 

we appeal to Lemma 9, which asserts the above inequality granted we show that uj^ A v e as 
N — > oo. To this end, observe that the differences j/j,fc+i — yj,k cannot be too large to often. 
Indeed, let Bn be the number of the differences Vj,k+i — Vj,k that are large than N^ 1 / 4 . Since 
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the sum of the heights of the rectangle Rj t k for a fixed index j is bounded by diam(supp(^ £ )) 
independently of j, it holds that 

+ c 2 Ar1/ ^ < diam(supp(f E ))M, 



jyi/4 * N i/2 

which immediately implies that Bn = O (A^ 3 / 4 ). Hence, if / is a function that is bounded by 1 in 
modulus and satisfies the Lipschitz condition with constant 1 on supp(z^ e ), then we can write 



M Mt 



M 



= E E / [/ - / (»&,*)] d ^ + / 

3 = \k=\ J jRC - j=l 

Observe that each of the last two terms on the right-hand side of the equality above is of order 
N^ 1 / 2 as |/| < 1, v e {R c ) < N' 1 / 2 , and M = pV 1 / 2 ] . Further, split the double sum into two: one 
over those rectangles that have heights at most TV -1 / 4 and the rest of them. Since the number 
of the "large" rectangles is B^ = O (iV 3 / 4 ) and |/| < 1, we get that this part of the sum is of 
order O (TV -1 / 4 ). On another hand, using the Lipschitz continuity of / on the first group of the 
rectangles and since the width of each rectangle is b^j/M, we deduce that 



f-f(<C)\du e = 0(N-^) 



and therefore the first part of the sum is also of order O (N 4 / 4 ). Hence, dist (v s , u)„N,e) = 
O (iV -1 / 4 ) and the claim follows. 

Finally, we define rj N ' e :— ' £ j by selecting any N point from the collection r\. As we are 

discarding at most N 1 / 2 points, uj v N, e still converges to v E in the weak* topology. Moreover, the 
behavior of the discrete energies also remains unaltered as the absolute value of the contribution 
of the removed points is of order O (N^ 1 / 2 log TV). 

Lemma 11. Let v e be defined by (13), r] N ' s be as in Lemma 10, and 
On.s ■ 



□ 



N,e 



< 



3VN 



x < z N 



N,e\ 
Z N - Vn I 



< 



3VN 



Then for alP tj € Oiv,ej it holds that 

dist (i/ e , lj v ) < O (AT 1 / 4 ) and 
for some constant c" that depends only on v e . 
Proof. It holds for every if £ On,e that 



| /* [Ur,] - I* [w^ 



< c"N- 1/2 log AT 



dist (a;,,, Wf.it,, 



sup 

/ 



fduJ v - / fduif, 



1 N 

3=1 



m | < 



3VN 



by Lipschitz continuity of /. Thus, the first claim follows from the triangle inequality. Furthermore, 
it holds that 



\I*[oj v ] - I*[uJ v N,e)\ < ^ E 

n^rn 



log 



1 



N,e 



N,e 



< 



1 



N,e< 



For brevity, we slightly abuse the notation and assume that r) 6 Ojv e stands for 7] = {zi, ■ ■ . , Zjv}, where 
(zi, . . . ,z N ) e Ojv, e . 
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Replacing the collection l 7 ?^' 6 } by {ffo.fc} constructed in Lemma 10, we shall only increase the 
sum on the right-hand side of the above chain of inequalities. As 

Mj M Mi M 

for a fixed point r}j,k, the desired result follows. □ 

Proof of Theorem 2. We start by proving the upper bound in (7). By Lemma 6, each extremal 
configuration {Ai,...,Ajv} realizing the supremum in (11) is, in fact, a configuration of Fekcte 
points for K since wk = 1 on K. That is, 

8n(w k ) = 5 N {K) := sup < ] [ \z n - z m \ : (zi, . . . , z N ) G K N \ . 

\m<n ) 



Hence, we get from (6) that 

at 



<4W (J/k +1 - N) ^ 
By (5) and since wk{z) ~ as \z\ — > oo, it holds that 



limsupA^log ( / w { " +1 ~ N}IS dA) < lira AT 1 log! / w 2 + c °dA ) = 0. 



j k J — Iy s> \ / K 

In other words, 

limsupAr- 2 logZ N , s . < [3 lim N~ 2 \ogd N (K) = ~I[w K ], 

AT-KX3 W->oo ^ 

where the last equality is a well known Fekete-Szego theorem [12, Thm. 5.5.2]. 

To prove the lower bound in (7), let K m be defined by (12). Further, let uj m ^ be the measure 
defined by (13) for v = ujK m and any e S (0, ^). Clearly, supp(w m . £ ) C K. Then we get from (6) 
that 



Z N ,s,p > I J] \Zn-Z m \ dA® N ( Zl ,...,Z N ) 

Jsupp(w m , e ) JV m<n 

= / Y\ \ Zn ~ Zm \ P CXP \ ~ E l°g a m, £ (^n) > dw®^(zi, . . . , 
m<n \ n— 1 J 



Hence, it follows from Jensen's inequality that 

logZjv, s ,/3 > -/(^Vlog-i - - + Vloga m , e (z n ) dw®^(2i,...,ZTv) 



7^e io stt4^+E ] 

J \ m<n \ Zn Zm < n=l 

N(N-1) , . /" 
-p ^ ^Kn>eJ --'V / loga m , £ da; n 



As a m)E < l/7re 2 , the integral J log a mj6 du; m . £ is finite and therefore 

liminf — - \obZm , « > — — lim Ifo;„, F l = — — I\u>k 1, 

where we used Lemma 8. Thus, the lower bound in (7) follows from Lemma 7. 
For brevity, put 

P N (v, e) := Prob {dist (v, uj v ) < e} , r\ = {qi, . . . , r/ N }. 
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Further, define 

O v>t := |(7?i,...,7?jv) G C N | dist (//,£«;,,) < e, 77 = {771, . . . , r/jv} 

Observe that 0„, e is an open subset of C N for any 1/ and e > 0. 

To prove the upper bound in (8), let a N ' e be a configuration that maximizes 

/ AT \ 



exp < -(/3s - 2 - cp) gxizn) \ [\z n - 



n — l ) m<n 



among all Appoint configurations a satisfying dist (i/,uj a ) < e. Such a configuration always exists 
since the function above decays as z -(fii s ~ N + 1 )~ 2 - c o) D y (5) m each coordinate. Therefore we are 
simply looking for a place where a continuous function reaches its maximum on some sufficiently 
large ball in C^. Then 



\ogP N {v,e) = log n N , Stf3 (z 1 ,...,z N )dA® N (z ll ...,z N )J 

< -\ogZ N , sJj -(0s -2- c )7V J g K doj aN , t - ^N 2 h[uo a N,] +N\og J 



g-(2+Co)ff_KT 



Let v e be a weak* limit point of {w q n, s } (again, as all the configurations a N,e belong to a ball of 
fixed radius in C N , the measures w a jv, e are compactly supported). Then, along the subsequence 
for which uj a N, t — > v e , it holds that 

limsupi\T a logP w (i/,e) < U I[u K ] - h[v<] 

V->oo * \ 

by Lemma 9, limit (7), and since J gK<iuj a N,e — > J g^diA Let now {^ ejl/ }^/ =1 , ejw — > as M — > 00, 
be an arbitrary sequence of weak* limit points constructed above. Since dist (y CM , v) < ejvf, it 
holds that v tM A v as M — > 00. Thus, the principle of descent and continuity of <?/<- yield that 
liminf m->-oo h[ v<iM ] > -^M- This includes the case supp(^) n K c 7^ and ^ = 0. In this situation 
supp(^ ejw ) C\ K c ^ for all large M and I [v] = Io[u eM ] = 00. Therefore, 

2, „ ,.. ' 



lim sup lim sup N logPjv(i/, e) < — I /[wjf] — I^] 

e->0 JV->oo 2 \ 

Observe also that the above considerations imply the existence of the full limit in (8) when 
I[v] = 00 or when £ = and supp(^) n K c ^ (clearly, the limit is —00). Thus, for the lower 
bound in (8), it is enough to consider measures with finite energy and, when £ = 0, only those that 
are supported on K. 

Assume first that I > 0. Let v be a Borel probability measure with finite energy and v t be the 
measure defined by (13) for some e > 0. Let further \r\ N ^ be a sequence of configurations con- 
structed in Lemma 10 and On,h be neighborhood constructed in Lemma 11. It follows immediately 
from Lemma 11, the triangle inequality, and Lemma 8 that On,e C 0„ ei£ C Ut 2e- Hence, 



P N (v,2e)>—^—[ cxp{-/3 S V 5if (z„)l TT 



\z n -z m fdA® N ( Zl ,...,z N ). 
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Then we deduce from Jensen's inequality that 

logP/vfy, 2e) > -logZAr, Si/ 3 + log|Oiv, ( 



> 



2 L, c [ N IM + 2sN J 9KdLu v — \o^\ — 

- log Z N . s . fs ~ | (n 2 I* [ Uv n..] + 2sN J g K duA + N log (4(c') 2 7r) 
-9AHog N - ^ (c"NlogN + 2sN fmax ^ g K dcJ v - J g K o\v, 



by Lemma 11 and since AA® N /|Ojv,e| is a probability measure on N . e , where |Ojv,e| = (4(c') 2 7r/9A^) 
is the volume of On,c- By compactness of On,* it holds that 

max J g K du)ri = f g K du) VN 

for some rj N — {f)i , ■ ■ ■ , r)j$} such that (t)i t ■ • • t Vn) ^ On,c- By Lemma 11, it holds that 

therefore 

max / gxdujr] — J gKdv e — > 



N 



(jj n — > v F as N — > oo and therefore 

'IN fc 



as TV —> oo. 



Hence, 



liminf TV" 2 log P N (v, 2e) > f ( J[wjf] - i>[i/ e ] 



by (7) and Lemma 10. The lower bound in (8) follows now from Lemma 8. 

Finally, assume that I = 0. Let v be a probability Borel measure with finite energy supported 
on K and {v m } be a sequence of measures granted by Lemma 7. Then (iVi)e is supported in K° 
for every e < l/2m and so are the measures lo v n,c and constructed in Lemmas 10 and 11 for 
(v m )e- Thus, the argument we used for the case £ > yields now that 

liminf liminf N~' 2 P N {v m ,£) > t;(i[wk] ~ hV'm] 

To show the lower bound in (8) it remains only to observe that Pn{v, e + e m ) > -P/v(*An, e), where 
e TO = dist (i/, v m j , and I[f m ] — > I[v] as to — >• oo by Lemma 7. □ 

Proof of Proposition 3. Let {^sl^o 1 ^ e the sequence of orthonormal polynomials with respect 
to the inner product 

(/,<?) = / /gcxp{-2 S5x }dA. 
Denote by x n ^ s the leading coefficient of ir n>s . Then it holds [5, Sec. 5.4] that 

N-l 

\ogZ N , s ,2 = log AH - 2 logxr„ iS . 

n=0 

It was shown in [16, Thm. 11 that K as described 



/n+ 1 ( n+l 



l + O 



n>s C p«+ i (iio y 7T 

Hence, it follows from (1) that 

log Zat !00 , 2 = -N(N + l)I[u K ] +N\ogTT + 0(l) 
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N-l 

log = log Z N>00>2 - lo s( s ~k)+Nlogs + 0(1). 

fe=0 

It follows from the Stirling's formula log rt! = (n + 1/2) logn — n + 0(1) that 

iV-l 

- ^ log(s - k)+ N log s = N+(s-N) log (l - s _1 .2V) + 0(logiV) 
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which yields (9). 
For e > 0, set 



□ 



(14) Un,e := I (zi, . . . , z N ) | Y[ \z n -z m \ ~[[w% 1 (z n ) >expj-(/ 
Lemma 12. For each N £ N and e > 0, it holds that 



N(N - 1) 



/ Sl NiSifl dA® N < exp l-/3{e + o(l)) ^ 1} } 

Jc N \u N . e I ^ J 



/C«\C/jv 

Proof. For (zi, . . . , zjv) € C w \ fjv e it holds that 



^iV,8,^(^l) - j z Jv) 



7V(iV - 1) 



N 

u 

n=l 



< exp(-/9(/[wjr] 

< exp{-/3(e + o(l))ffi-^)n^ +C0 (^) 



Ms-N+l) 
J K 



( z n) 



where we used (7) and (5) for the second inequality. Now, the conclusion of the lemma follows 
from the fact that 

N / „ \ N 



I IT w 2 + C0 (z n )dA» N <( f w 2 K +Co dA) = exp {O(N)}. 

Jc N \u N , c n=1 \Jc J 



□ 



For / G C fc (C n ), set 



(15) 



/n(zi, ...,z n ) 



(N — n)\ 



where the sum is taken over all distinct permutations a of size n for {1, . . . , N}. 

Lemma 13. Let f G Cf,(C") and z = {zi, . . . , z^r}. Then there exists a finite constant c(n) such 
that 



/n(zi, ...,z n )- / fdu 



< max |/| 



c(n) 
TV 



Proof. For ?i = 1 it simply holds that fxizi, ■ ■ ■ , Zn) = J fdui z . When n = 2, it is true that 



/at(zi, . . . ,zjv) - / fduj 



1 



TV- 1 



/dw 



®2 



N - 1 



/(u,ii)du; z (ii) 



< 



2maxo | /| 
AT- 1 



More generally, it holds that /jv(zi, . . . , Zn) is equal to 
(N -n)\ 



N n J fdujf " - iV"- 1 £ / /d^f"- 1 + ■ ■ • + (-l)"" 1 ^ J f(u, . . . , U )d^( t 
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where the first sum is taken over all possible combinations of two coordinates being equal, the next 
sum is taken over over all possible combinations of three coordinates being equal, etc. Since the 
number of terms in each sum depends on n but is independent of N, the conclusion of the lemma 
follows. □ 

Denote by C c (C n ) the collection of continuous functions on C™ with compact support. 
Lemma 14. For any f € C C (C"), it holds that 



lim lim 



f^O N- 



sup 

(z 1 ,...,z N )eu N 



,z N )- I fdujf 1 



= 0. 



Proof. For fixed e, let | ^Af ' £ , . . . , A^' £ ^ | . A C M, be a maximizing sequence for the first limit 

in question. Since / has compact support, we can assume that |A^' £ | < R for some R = R(f) 
large enough. Set 

Let v e be a weak* limit point of {lj x n.e}. Clearly, i> £ is also supported in the disk of radius R. 
Then it follows from Lemma 13, the choice of j^A^' 5 j that 



lim sup sup 

IV— >0O (zi,...,2jv)GfA 



f N (zi,...,z N ) - / fduj 



fdcj 



K 



(observe that — > vf n since lu x n, s A v £ ). Further, let v be a weak* limit point of {i> £ } such 
that 



lim sup lim sup 

N ^°°(z u ...,z N )£U N , e 



f N (zi, - ..,«jv) 



fdv 



fduj 



SO;) 
K 



The proof of the lemma will be completed if we show that v = ujk- To this end, recall that 
(^X^' £ , . . • , A^' 6 ^ <E Un, s and therefore 

I*[cl) a jv,e] + 2 J g K dui X N, c < I[wk] + £• 

Then it follows from Proposition 1, Lemma 9, and the inequality above that 

I[o>k] <h\v e ] </[wjc]+e. 
Applying principle of descent once more, we get that 

I[w K ] < h[v] < liminf h[i/ e ] < I[u K ]- 
The desired conclusion now follows from Proposition 1. □ 

Pivof of Theorem 4- Fix / £ C c (C n ). That is, we assume that / has compact support. Define /at 
by (15). Then 



c N 



f N n N , St pdA® N = 



(N-n)\ 



/ / (Mi)' ■ • ■ ' M™)) ^n.s,/3(zi, z N )dA® N (z ll ...,z N ) 
„ Jc N 



C" 
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as £In,s,/3 is symmetric with respect to the permutations of z\, . . . , Zn. Moreover, since £In,s,p is a 
probability density function by its very definition, it holds that 

f fn^gdA®"- [ /duf"= / ( f N - [ fduj® n ) n N , s<p dA® N . 

Hence, by Lemma 12 we only need to show that 

lim / (fN - [ fdujf 1 ) n N . s . p dA® N = 



for any e > as sup c « |/at| = sup c „ |/|. The desired conclusion now follows from Lemma 14 since 



In - fdu) 



A" 



f N (z 1 , ...,z N )- I /dwf 



Cl N>St pdA® N < sup 

(z 1 ,...,z N )GU N ,. 

where we once more used the fact that Qn,s,p is positive and has unit integral over C . 

Now, let D be a large enough ball in C™ to contain K n and / be a function supported in D 
satisfying < / < 1 and such that / = 1 on K n . Then 

lim sup / n% l \ g dA® n < 1 - lim / fSl™ g dA® n = 1 - / fdujf 1 = 0. 

Since any / g C&(C") can be written as a sum f c + (/ — / c ), where f c € C C (C") and / = f c in D, 
the general claim follows. □ 
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